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The pair-potential and current density around a single vortex of the two-dimensional chiral 
p-wave superconductor with d = z {px ± ipy) are determined self-consistently within the qua- 
siclassical theory of superconductivity. Shrinking of the vortex core at low temperatures are 
considered numerically and analytically. Temperature-dependences of the spatial variation of 
pair-potential and circular current around the core and density of states at zero energy are the 
same as those in the isotropic s-wave case. When the senses of vorticity and chirality are oppo- 
site, however, we find two novel results; 1) the scattering rate due to non-magnetic impurities 
is considerably suppressed, compared to that in the s-wave vortex. From this observation, we 
expect that the chiral p-wave superconductors provide the best chance to observe the shrinking 
of the vortex ("Kramer-Pesch effect") experimentally. 2) The pair-potential of chiral p-wave 
superconductors inside vortex core recovers a combined time-reversal-Gauge symmetry, although 
this symmetry is broken in the region far from the vortex core. This local recovery of symmetry 
leads to the suppression of the impurity effect inside vortex core. 
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§1. Introduction 

The superconducting state in Sr2Ru04^' is expected to 
be unconventional; a Knight shift measurement in pla- 
nar magnetic field suggests odd-parity pairing with d 
vector perpendicular to the Ru02 plane. ^' Muon spin 
relaxation experiment suggests spontaneous magnetiza- 
tion;^' this means that time-reversal symmetry is bro- 
ken spontaneously in the superconducting state. The 
simplest pairing symmetry which is consistent with the 
above two experimental results and the tetragonal crys- 
tal symmetry is given by the odd-parity pairing with 
d = z{px iipj,).^' This is commonly called the "chi- 
ral p-wave superconducting state" because the Cooper 
pair has a definite angular momentum L^ = ±1 in the 
internal motion. A singly quantized vortex is, on the 
other hand, a configuration of the pair-potential where 
the center-of-mass motion of Cooper pair has an angular 
momentum. The interplay between vorticity and chiral- 
ity would bring about intriguing phenomena in vortex 
physics of Sr2Ru04. Indeed, several authors^"^-' have ar- 
gued that the impurity effects in vortex core of chiral 
p-wave superconductor is considerably suppressed, com- 
pared to that in s-wave vortex core. The fact that the 
superclean regime is realized within the core of chiral p- 
wave vortex has several implications. First, the fiux flow 
conductivity is expected to be enhanced at low temper- 
atures in the chiral p-wave vortex states.®' Second, a 
large Hall effect is also expected. Third, the shrinkage 
of vortex cores at low temperatures, which is referred 
to as "the Kramer-Pesch (KP) effect" ,^' is free from the 
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smearing due to impurities. The vortex shrinking effect 
in chiral p-wave superconductors is the main issue in this 
paper. 

The KP effect, originally predicted for s-wave vor- 
tices, was studied numerically by Gygi and Schhiter,^*^' 
who gave a self-consistent solution to the Bogoliubov- 
dcGennes (BdG) equation. More general cases were dis- 
cussed by Volovik.^^' Vortex shrinking in d-wave case 
was confirmed numerically by Ichioka et al}"^' Exper- 
imental studies on the vortex shrinking in NbSe2 and 
cuprates have been reviewed recently by Sonier et al}^> 
For the chiral p-wave case, the electronic structure and 
charging effect were studied^^^^^' within the framework 
of the Bogoliubov-deGennes (BdG) equation. While the 
analysis in the quantum limit k^S, ~ 1 brings us the un- 
derstanding of physics at extremely low temperatures, 
where the quantum effects becomes important, ^^"^"■' the 
quasiclassical condition fcpC ^ 1 is satisfied in Sr2Ru04 
with a ^0 (the coherence length at low temperature) 
« 660A.^^' In this paper, we will study the structure 
of chiral p-wave vortices by adopting a more appropriate 
framework; the quasiclassical theory of superconductiv- 
ity. ^^"^^' The quasiclassical scheme has previously been 
applied to the study of the pair-potential and the local 
density of states near domain wall.^^' 

In the next section, we summarize the formulation of 
the quasiclassical theory of superconductivity. In § 3, 
we calculate the pair-potential analytically, only consid- 
ering contributions to the gap equation from the chiral 
branch. In § 4, we present our numerical results for the 
pair-potential and establish the validity of our analytical 
findings of § 3. In § 5, we discuss aspects different from 
s-wave vortex cores. In S 6, we summarize our results. 
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§2. Formulation 

2. 1 Quasiclassical theory of superconductivity 

We consider a two-dimensional superconductor with 
the circular Fermi surface. For simplicity, we consider 
the type II limit (k = X/S, -^ cxd) in this paper except 
in subsection 3.4. In Sr2Ru04, the coherence length ^o 
at zero temperature is expected to be « 660A^^) and 
much larger than the Fermi wave length, which is of the 
order of several angstroms. Since the atomic-scale struc- 
tures are irrelevant in determining the pair-potential, we 
can factor out the rapidly oscillating terms which enter 
the wavcfunctions of the BdG equation. The resulting 
equation (Andrecv cquation)^^-* is equivalent to the qua- 
siclassical equation of superconductivity, ^^^^'''•' which we 
now explain briefly. The quasiclassical theory of super- 
conductivity is described in the equilibrium case in terms 
of the quasiclassical green function 



giz,r,p) ^g{z,r,a) 



(2.1) 



9 f 
. -f -9 

which is a 2 X 2 matrix in particle-hole space and 
is a function of (complex) frequency z, the direction 
p — (cos a, sin a) of momentum p = pfP, a point 
r = r(cos (j), sin 0) in real space. The equation of mo- 
tion for g is given by 



-iv{p) -Vg = ef3-A(r,p) 
supplemented by the normalization condition 

g = -TT 1, 

where 1 is the 2x2 unit matrix. 
Here fs is a Pauli matrix, 

1 
-1 



and A is given by 



A(r,p) 
-A*(r,p) 



(2.2) 
(2.3) 

(2.4) 
(2.5) 



where A{r,p) is the pair-potential. 

To determine the pair-potential self-consistcntly, we 
have to consider the case where z is a Matsubara fre- 
quencies iujjn = IttT (2m -I- 1). Setting z = e + id with 
real e, on the other hand, we can obtain the density of 
states from —Inig{e + i5). 

The self-consistent equation for the pair-potential 



A(r, a) = A+ (r) e'" + A" (r) e" 



(2.6) 



is given by 
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A±(r)^Ty, J2 / 
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da 

2^ 



e^'"f{icj^,r,a). (2.7) 



Here the strength Vp of the attractive interaction is given 

by 



Wc/(27rT) 



V-'^ In {T/T,)+ J2 



1 



m+1/2' 



The current density j{r) is given by 

Ef da ^ 



(2.8) 
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where A'o denotes the normal state density of states and 
e(< 0) charge unit. The set of equations (2.1)-(2.7) pos- 
sess two continuous symmetries, the rotational symme- 
try around z-axis with the generator L^ and U(l) Gauge 
symmetry with /. The actions of these generators are 
written, respectively, by 



and 



L, = -i {8/ del) + 8/ da) 
/A = A, /A* = -A*, Ig^O 

if = f, if = -I 



(2.9) 



(2.10) 



We note that the expression (2.9) consists of the center- 
of-mass rotation {—\d/8<j)) and the internal rotations of 
the Cooper pairs {—id /8a). 

We consider two configurations of a single vortex with 
the following asymptotic forms for the pair-potential: 



iW+a 



Case I 
Case II 



(2.11) 



A(r. ,._. , ^^^^ 

and for the green function: 

^l X T^ ( -z, A{r,a) 

^^'' "■' "^ " VA(r,a)2-z2 1, -A*(r, a), 

(2.12) 
In (2.11), Aoo denotes the amplitude of the pair-potential 
in the bulk and is given by 

,2 \-l/2 



\u>m\<UJa 



The asymptotic forms (2.11) and (2.12) determine the 
symmetries of the two vortices; each symmetry is de- 
scribed by the modified generator Q = Lz ~ ml with 
m = for case I and m — 2 for case II. 

By solving QA{r, a) = 0, the general expression for 
pair-potential in case I is obtained as 

A(r,a) = A+(r)e'(-'^+")+Ao(r)e'('^-") (2.13) 

and that in case II is given by 

A(r,a) = A+(r)e'('^+") + Ao (r)e'(3<^-"\ (2.14) 

Here Aq (r) denote functions of r — \r\. In the next 
section, we will numerically determine the self-consistent 
solution for Aq (r). 

2.2 Riccati Formalism 

In order to determine the pair-potential, we could solve 
the quasiclassical equation (2.2) and the gap equation 
(2.7). Instead of directly attacking eq. (2.2), we re- 
formulate the original quasiclassical equation (2.2) into 
Riccati-type equations for auxiliary functions. ^*"^'^'' Let 
7(2:, r, a) and 7^(z, r, a) be the solutions of the following 
Riccati differential equations: 



IV ■ V7 



A * 2 

-A 7 



2e7- A 



(2.15a) 
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iv ■ V-f^ = -A7^2 + 2€-f^ - A*. (2.15b) 

The solution g of (2.2) can then be given by 



9 = 



1 — 77^ 27 



l+77t I 27t -1 + 77+ 



(2.16) 



In this parameterization, the norniahzation condition 
(2.3) is automatically satisfied. 

For a given a, the framing of space 

r ^ sp + bz X p 

is naturally introduced. Since the Fermi surface is circu- 
lar symmetric, v{p) ■ V can be rewritten as 

vp ■ "V = vd/ds. 

Thus, eq. (2.15) reduces to a one-dimensional problem 
on a straight line (quasiclassical trajectory) with a con- 
stant impact parameter b = rsin((j) — ct). Equations 
(2.15) are to be solved, under the initial conditions: 

7 = 1 I, for s — — Sc (2.17a) 



for s = Sc, (2.17b) 



-e-i(5-i^|Ar-e' 
A* 



e + iS + h \Ar -e 



where Sc is a positive real cut-off. Numerical calculation 
of 7 for Imz > is stable if (2.15 a) is integrated start- 
ing with (2.17 a) towards increasing s. The calculation 
of 7I for Imz > 0, on the other hand, is numerically 
stable if eq. (2.15 b) with the initial condition (2.17 b) 
is solved with decreasing s.^^^ Obviously eq. (2.7) con- 
tains /(iciJm, r. Of) with negative Wm- Such contributions 
can, however, be eliminated easily with the use of the 
symmetry relation /(io^m, r, a) = —f{—iuJm,'r',a + TT).'^^' 
Numerical integrations of (2.15) are performed by the 
adaptive stepsize control Runge-Kutta method. ■^^^ The 
cut-off frequency luc is taken as lOAoo, where Aoo is the 
modulus of the pair-potential away from vortex center 
at zero temperature. The cut-off length Sc is taken as 
3'Coi 6^0; 12^0 (Co = v/Aoo)- The choice of cut-off de- 
pends on temperatures. 

§3. Analytical Results 

3. 1 Chiral branch and Kramer-Pesch effect 

In both cases (I and II), each vortex has a zero-energy 
bound state for trajectory with 6 = 0, because the order 
parameter is real up to an overall phase and changes its 
sign when going from s ~ —00 to s = 00."^^' On a tra- 
jectory with nonzero impact parameter, there is at least 
one Andreev bound state with nonzero energy. The low- 
est bound states for each impact parameter form a sin- 
gle branch, which corresponds to the Caroli-deGennes- 
Matricon mode.'^^^ Roughly speaking, this branch has 
energy linear in the impact parameter b ("chiral branch 
" in this sense) and merges into the bulk continuum spec- 
trum with & ^ Co- Therefore, the dispersion can be es- 
timated as -E '-~^ Aoob/^Q. This chiral branch plays an 
important role in the self-consistent calculation of the 
pair-potential, as explained below. 

Overall features of the pair-potential around a single 



vortex at a wide range of temperatures lower than Tc can 
be described by the GL theory; Although the GL theory 
is valid only near Tc in principle, the length scale which 
the bulk quasiparticle introduces into the gap equation 
(2.7) is of the order of Co even away from Tc. However, 
at temperature T(<C Tc), the chiral branch with energy 
comparable to T gives a large contribution to the pair- 
potential in the gap equation (2.7) at the radius r ^ 
^1 = T^o/^00 from the vortex center. Thus, the chiral 
branch introduces another length which scales the spatial 
variation of pair-potential in the vicinity of vortex center. 
The emergence of ^1 at low temperatures can be regarded 
as a result of the shrinking of vortex core (KP effect^)). 
For the chiral branch with |z| ^ Aqo, the green func- 
tion g can be obtained analytically.^' ^'"^^^ The contri- 
bution in the gap equation (2.7) from the chiral branch 
can, therefore, be calculated analytically. Further we as- 
sume that near the vortex center, contributions in (2.7) 
to the pair-potential from other excitations are negligi- 
ble. Within this assumption, the self-consistent deter- 
mination of the pair-potential becomes feasible. This 
assumption will be justified by the numerical calculation 
in §4. 

3.2 Case I: A(r, a) -> Aooe'^"^"") 

The analytic expression for the green function of the 
chiral branch in the s-wave case was first obtained in 
ref. (9). We, however, follow the method presented in 
ref. (33), which is based on the Riccati formalism. The 
outline of the calculation of g is as follows: First, calcu- 
late 7 and 7+ for z = and b = Q. This nonlinear dif- 
ferential equation happens to be elementarily integrable. 
Second, calculate 7 and 7' up to linear order in z and 
b. The problem reduces to a first-order linear differential 
equation, which can be easily solved. Third, substitute 
the resulting 7 and 7+ into (2.16). The zero-th order 
with respect to z and b is cancelled in the denominator 
in (2.16). Therefore, the first order term should be kept 
there. In the matrix part in (2.16), on the other hand, 
zero-th order term with z and b survives and hence the 
first order contribution can be neglected. One can refer 
to refs. (8) and (33), for the details of the method. 

The resulting expression for g{z, r, a) is given by 



TTu exp [— u_ (s)] M_ 
2C-{z-E^{b)) ' 



(3.1) 



with 



M_ = 



1, 



Here exp [— m_(s)] with 

u-(s) = -/'''ds'{A„-(s') + A+(s')} 
^ Jo 

is essentially the spectral weight (up to a constant factor 
) of the Andreev bound state for the quasiclassical tra- 
jectory specified by b and a. The function exp [— 'u_(s)] 
takes the value of the order of unity for s <C Co a-nd much 
smaller than unity for s ^ ^o- The constant C_ in the 
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denominator of (3.1) is defined by 

/•OC 

C_ = / dscxp[-u_(s)]. (3.2) 

Jo 

The symbol E^{b) denotes the energy of the Andreev 
bound state, which is given by 

h_ f°° ds 
IT- 



E-ib) = — I -exp[-u^is)]{A-is')^A+{s')}. 
Jo * , , 

(3.3) 

Substituting the expression for f(z,r,a) in (3.1) into 

(2.7), we obtain 

T7:vVph{E'-r/{2T)) 



Ao*(r) 



4C- 



(3.4) 



with 

In{R) 



"^ da' 



sinna' tanh (_RsinQ;') n=l,3, • 



Here E'_ is defined by dE^{b)/db. For r < ^i_ 
2T/E'_, the expression (3.4) reduces to 



(3.5) 



Ao^(r) 



8C- Ui- 



(3.6) 



From eq. (3.6), we can see that both components Aq and 
A^ vanish linearly with respect to the distance r from 
the vortex center. Further, both components have the 
same amplitude. The relative phase between them is tt. 
The initial slope of Aq (r) is shown in Fig. 1. 

The length ^i_ can be estimated in the following way. 
In (3.2) and (3.3), exp [— ■u_(s)] becomes unity for s <C Co 
and much smaller than unity for s 3> ■fo- Thus C_ in 
(3.2) is of the order of ^g and then E'_ is reduced to 

From the above expression, we obtain 

ei_/eo = 2T/{CoEL) ^ (r/A^)/ln(eo/Ci-) , 

which gives 6-/Co - (T/Aoo)/ln(Aoo/r). According 
to this argument, the slope of the pair-potential around 
the vortex center is found to diverge in the limit of zero 
temperature. 

The current density coming from the chiral branch 
(3.1) is given by _7_(r) = j^{r)(j) with 



J- (r) 



2cr AC- V6^ 



-Jc(0) 



r 



(3.7) 



Here jc = 2|e|AfovAoo denotes the critical current in a 
homogeneous state at zero temperature. The contribu- 
tion to the current density from the extended quasipar- 
ticles are expected to be of the order of jdr/^o)- At 
temperatures much lower than Tc, where Ci_ <^ £,o, the 
distribution of the current density becomes more concen- 
trated around the vortex center. 

The shrinkage of the vortex core also affects the den- 
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Fig. 1. Analytical results (eq. (3.6) in the text) for case I where 
total angular momentum is zero. The initial slopes of the pair- 
potentials (a) Ar (r) for the dominant component with negative 
chirality and (b) the induced component Aq (r) with positive 
chirality are shown as functions of the distance r from the vortex 
center, for temperatures T/Tc = 0.02 (solid line), 0.05 (dashed 
line) and 0.1 (dotted-dashed line). The vertical and horizontal 
axes are scaled, respectively, by the modulus of the bulk pair- 
potential AcxD and the coherence length ^o a-t zero temperature. 



sity of states (integrated over a vortex core) : 
N{e = 0) = -— f — f drlmg{z = + [6). (3.8) 

TT 7o I" J 

Substituting (3.1) into (3.8), we obtain 
N{e = 0) « TVoCoV In (^0/6-) - A^oCoV In (Tc/T) . (3.9) 

This logarithmic factor (1/ In (Tc/T)) can be seen, in 
principle, in the specific heat at low temperatures. 

3.3 Case II: A(r, a) -^ Aooe'(^+") 

Following the method in refs. (33), (8), the expression 
for g in the present case is given by 



TTV exp [— U-(- (s)] M_|_ (a) 



2C+{z-E+{b)) 



(3.10) 



with 



M+{a) = 



o2iQ 



-2iQ 
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Here the spectral weight exp [— m+(s)] is given by 

"+(^) = -/'d.'{A+(.') + A„- (.')}. 
" Jo 

The constant C+ is defined by 



/ dsexp [— u+(s)] . 



The roles of m+(s) and (7+ are similar to those of m-(s) 
and C_ in the previous subsection. 

The dispersion E+{h) of the Andreev bound state is 
given by 

7 /*C>0 1 

E+{h)^— -expHi+(,s)]{A+(,s') + 3A-(.')}. 

^+ Jo s , ^ 

(3.11) 

From (3.10) and (2.7), we obtain 



Ao+(r) 



TTUIC 



4C4 



"h 






r 



(3.12a) 



-— . (3.12b) 



Here ^i+ is 2T/E'^ with E'_^ = d^+(fe)/d6. For r < 
^1+ = 2T/i?^, the expression (3.12) reduces to 



A+(t 



A-(r) 



8C+ 



TTuVn 



r 



(3.13a) 



96C+ va 



-— . (3.13b) 



The result (3.13) is shown in Fig. 2. The relative phase 
between the two components in this case is zero. If we 
scale the r-dependence in (3.13) by a+j the ratio of the 
prefactor of each component then takes a universal value 
{— 12). From (2.8) and (3.10), the current density is 
given by J_|_(^) = j+{r)(f> with 



J+W 



-Kev'^Noh (7'/a+) ttcv^Nq / r 



2C: 



4C4 



a 



(3.14) 



By an argument similar to that in the previous subsec- 
tion, we find that C+ ^ ^Oi 



and 



i?;^(Aoo/a)in(a/a+) 



a+/a-^(r/A^)/ln(A^/T). 



Temperature dependence of the length ^i+ leads to the 
shrinkage of the distributions of the pair-potential and 
circular current around vortex core at T much lower than 
Tc. 

3-4 Modification due to finite k 

In this subsection, we briefly discuss on the case where 
the GL parameter k is finite; most results in the previous 
subsections hold with the following modifications. 

In section 2, e in eqs. (2.2), (2.15) and (2.17) is re- 
placed by e = e + evp ■ A{r)/c in the case of finite k. 
The symbol c denotes the velocity of light. Owing to 
this change, E^{b) in (3.3), (3.4) and (3.6) is replaced 
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Fig. 2. Analytical results (eq. (3.13) in the text) for case II where 
total angular momentum is two. The initial slopes of the pair- 
potentials (a) A|f (r) for the dominant component with positive 
chirality and (b) the induced component A^(r) with negative 
chirality are shown as functions of the distance r from the vortex 
center, for temperatures T/Tc = 0.02 (solid line), 0.05 (dashed 
line) and 0.1 (dotted-dashed line). The vertical and horizontal 
axes are scaled, respectively, by the modulus of the bulk pair- 
potential AcxD and the coherence length ^o ^^ zero temperature. 



hy E^{b) =E'_b with 
E'_=E'_ 



cC- 



dr 

— A(r)exp[-u_(r)]. (3.15) 



Here A{r) is introduced through A{r) = A{r) (z x f). 
The length scale ^i_ in (3.6) is replaced by 2T/E'_, ac- 
cordingly. Similarly, E^(b) in subsection 3.3 is replaced 
hy E+{b) =E'_^b, with 



E'=E' 



cC^ 



dr 



A(r)exp[-M+(r)]. (3.16) 



The length ^i+ is replaced by 2T/E'^. 

The current density coming from the chiral branch 
(3.1) is given by j{r) ~ j{r)<p with 



Jir) 



2C_ 



nev'^No 
4C- 



^)- (3.17) 



As in the previous subsection, the density of states (inte- 
grated over a vortex core) at zero energy (3.8) is obtained 
asNoCi/ln{T,/T). 
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§4. Numerical Results 

4.1 Case I: A{r,a) -^ Aoog'^"^"") 

Figure 3 shows the numerically determined self- 
consistent pair-potentials as a function of r/^o , at tem- 
peratures T/Tc = 0.02 (solid line), 0.3 (dashed line) and 
0.6 (dotted-dashed line). The relative phase between 
the dominant (A(^(r)) and the induced (Aq (r)) com- 
ponents is TT. Further, the spatial variations near the 
vortex center becomes steeper with decreasing tempera- 
ture and hence the shrinking of the vortex core at low 
temperatures is confirmed numerically in the chiral p- 
wave superconductors. These two points have been ex- 
pected from the analytical results in the previous sec- 
tion. We compare the numerical results at T/Tc = 0.02 
with the analytical result eq. (3.6) in fig. 4. For both 
components of the pair-potential, analytical results ex- 
plain the numerical results within the accuracy of 12% 
at r/^o = 0.0025. Figure 5 shows numerical results on 
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Fig. 3. Numerically determined self-consistent pair-potentials for 
case I where total angular momentum is zero, (a) Aq (r) of the 
dominant component with negative chirality and (b) the induced 
component Aq (r) with positive chirality are shown as functions 
of the distance r from the vortex center, for temperatures T/Tc = 
0.02 (solid line), 0.3 (dashed line) and 0.6 (dotted-dashed line). 
The vertical and horizontal axes are scaled, respectively, by the 
modulus of the bulk pair-potential A 00 and the coherence length 
^0 at zero temperature. 



the circular current density surrounding the vortex core 
for temperatures T/Tc = 0.02 (solid line), 0.3 (dashed 
line) and 0.6 (dotted-dashed line). The current density 
and the distance from the vortex center are, respectively. 
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Fig. 4. Comparison between analytical (solid line) and numer- 
ical (dashed line) results on the pair-potential for case I at 
T/Tc = 0.02. (a) A|7(r) of the dominant component with neg- 
ative chirality and (b) the induced component Aq (r) with pos- 
itive chirality are shown as functions of the distance r from the 
vortex center. The vertical and horizontal axes are scaled, re- 
spectively, by the modulus of the bulk pair-potential Aoo and 
the coherence length ^0 at zero temperature. In both figures, 
the relative deviations of analytical and numerical results are 
0.12 at r/^o =0.0025. 



scaled by the critical current density jc — 27Vo|e|i;Tc in 
the bulk and the coherence length ^0 at zero temperature. 
We note that the (fi component is always negative, as ex- 
pected from the analytical calculations. We can see in 
the figure that the distribution of the circular current be- 
comes more concentrated with decreasing temperature. 
This feature is consistent with the results on the pair- 
potential, which shows the shrinkage of the core at low 
temperatures. Figure 6 blows up the vicinity of the core 
center at T/Tc = 0.02. The numerical result is shown by 
the dashed line. The analytical result (solid line) is dis- 
played for comparison. Obviously, the numerical result 
approaches the analytical result asymptotically near the 
vortex core. 

4.2 Case H: A(r, a) -^ Aooe'('^+") 

Figure 7 shows numerically the determined self- 
consistent pair-potentials at temperatures T/Tc = 0.02 
(solid line), 0.3 (dashed line) and 0.6 (dotted-dashed 
line). Once again, we see that the core shrinking oc- 
curs even in this case. The relative phase between the 
two components of the pair-potential is in this case. 
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Fig. 5. Numerically determined current density is shown as a 
function of the distance r from the vortex center for case I at 
temperatures T/Tc = 0.02 (solid line), 0.3 (dashed line) and 0.6 
(dotted-dashed line). The vertical and horizontal axes are scaled, 
respectively, by the critical current density jc = 2No \ e \ vTc in the 
bulk and the coherence length fo at zero temperature. 
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Fig. 6. Comparison betvifcen analytical (solid line) and numerical 
(dashed line) results on the current density for case I at T/Tc = 
0.02. The vertical and horizontal axes are scaled, respectively, 
by the critical current density jc = 2No\e\vTc in the bulk and 
the coherence length go- 



These two points have been anticipated from the ana- 
lytical calculation. For the induced component AQ"(r), 
the spatial variation is expected to be proportional to 
(r/^o)^ near the core at low temperatures. Correspond- 
ingly, figure 7(b) shows the positive curvature near the 
core. Precisely speaking, we have confirmed that A^(r) 
at T/Tc ~ 0.1 is proportional to (r/^o)"^ in the range 
r/^o < 0.01. At the lowest temperature T/T^ = 0.02, 
however, the cubic dependence on r/S,o cannot be con- 
firmed at the distance larger than r/^o = 0.002 within 
our numerical accuracy. Thus, a comparison with ana- 
lytical result is made only for the dominant component. 
Figure 8 shows the analytical (solid line) and the numer- 
ical (dashed line) results on the pair-potential for the 
dominant component at T/Tc ~ 0.02. We find good 
agreement between the two results for distances smaller 
than r/^o = 0.0025. 

Now we turn to the circular current. In fig. 9, the 




Fig. 7. Numerically determined self-consistent pair-potentials for 
case II where total angular momentum is two. The pair- 
potentials (a) Aq (r) of the dominant component with positive 
chirality and (b) the induced component A^(r) with negative 
chirality are shown as functions of the distance r from the vor- 
tex center, for temperatures T/Tc = 0.02 (solid line), 0.3 (dashed 
line) and 0.6 (dotted-dashed line). The vertical and horizontal 
axes are scaled, respectively, by the modulus of the bulk pair- 
potential AcxD and the coherence length go at zero temperature. 
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Fig. 8. Comparison between analytical (solid line) and numerical 
(dashed line) results on the pair-potential Ag (r) for the dom- 
inant component with positive chirality for case II at T/Tc = 
0.02. The vertical and horizontal axes are scaled, respectively, 
by the modulus of the bulk pair-potential A^o and the coherence 
length go at zero temperature. 
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spatial distribution of the circular current is shown for 
T/Tc = 0.02 (solid Hne), 0.3 (dashed hue) and 0.6 
(dotted-dashed line). The unit of the current density is 
again jc — 2Nf)\e\vTc. Overall features are similar with 
fig. 5; narrowing of the distribution at low temperatures 
are equally observed in both cases. The magnitude of 
current density near the core in the present case is, how- 
ever, about 1.5 times larger than in case I. This means 
that the cost of the kinetic energy of the condensate due 
to the introduction of vortex is larger in case II than that 
in case I. This energy difference favors case I against case 
II as a stable state in the presence of magnetic field. This 
observation is consistent with the Ginzburg-Landau the- 
ory for chiral p-wave superconductors. ^^^ 

1. 




Fig. 9. Numerically determined eurrent density is shown as a 
function of the distance r from the vortex center for case II at 
temperatures T/Tc = 0.02 (solid line), 0.3 (dashed line) and 0.6 
(dotted-dashed line). The vertical and horizontal axes are scaled, 
respectively, by the critical current density jc = 27Vo |e|?;Tc in the 
bulk and the coherence length ^o a-t zero temperature. 



In fig. 10, we compare the numerical result (dashed 
line) with the analytical one (solid line) at T/Tc = 0.02. 
The numerical result approaches asymptotically and col- 
lapses on the analytical result at r/^o ^ 0.001. 

§5. Discussion 

Most aspects of the shrinkage of vortex core in chiral 
superconductors are found to be similar to those in the 
isotropic s-wave superconductors. However, there are 
two noticeable differences between the two cases, which 
is discussed below. 

First, we discuss the stability of the vortex shrinkage 
against the impurities. In refs. 5-8, it was shown that 
the impurity effects within vortex core in chiral super- 
conductors arc quite different from those in the isotropic 
superconductors. In the isotropic s-wave superconduc- 
tors, the impurity scattering rate T{e,T) for the chiral 
branch with energy e at temperature T was estimated to 
be r„ln(rc/e)/ln(Tc/T)38.39) in the Born limit, where 
Fn is the scattering rate in the normal state. The scat- 
tering rates F(e, T) for chiral superconductors can be es- 
timated by following the method used in ref. 8. In case 
I, F(e,r)/F„ turns out to be at most 0(e/Tc), while 
F(e,T)/F„ - l/ln(Tc/T) in case II. We note that the 
results in ref. 8 on F still hold in the present case al- 
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Fig. 10. Comparison between analytical (solid line) and numer- 
ical (dashed line) results on the current density for case II at 
T/Tc = 0.02. The vertical and horizontal axes are scaled, re- 
spectively, by the critical current density jc = 'iNo\e\vTc in the 
bulk and the coherence length ^o at zero temperature. 



though the induced (minor ) component has been ne- 
glected in ref. 8. In the presence of impurities, the length 
^i± should be replaced by 

a±=Max[T,F]/i?;(T). 

Thus ^i± becomes temperature-independent at T{< F) 
and hence the vortex ceases to shrink. Within the loga- 
rithmic accuracy, F ~ F^ in the isotropic s-wave and chi- 
ral p-wave (Case II) superconductors. In the two cases, 
therefore, the temperature range where the KP effect oc- 
curs is given by Fn < T <C Tc- In the chiral p-wave (Case 
I) superconductors, on the other hand, F is smaller than 
T under the condition T^/T^ <C 1 (clean superconduc- 
tors). Thus, in case I, the presence of impurity imposes 
no restrictions on the temperature range where the KP 
effect occurs. The importance of this fact will be under- 
stood by recalling the current experimental situations 
on the KP effect. ^"^^ The /i SR experiments on NbSc2 
with Tc w 7. IK show that the vortex core size ceases to 
shrink at T « IK and the temperature dependence of 
the core size at T > IK is much weaker than theoret- 
ical prediction.'"^) The saturation T « IK is presum- 
ably due to the impurity scattering, which smears out 
the structure of the energy spectrum of the bound states 
inside core. Takagi et al'^'^^ measured the flux flow con- 
ductivity of NbSe2 with high purity {RRR = 100) and 
did not observe the logarithmic temperature dependence, 
which could be regarded as an evidence for the KP effect. 
Therefore, if we try to observe the KP effect in NbSe2, 
we have to prepare extremely pure samples with RRR 
^ 100. Instead, chiral p-wave superconductors give us 
the opportunity to see the KP effect, by suppressing the 
impurity scattering rate in vortex core. 

Second, we make a remark on the local recovery of 
the symmetry in the pair-potential in case I (L^ = 0). 
In this case, the pair-potential (2.6) with (3.6) is in- 
variant under the combined action of time-reversal op- 
eration (A(r,a) — > — A*(r,Q; -I- tt)) (minus sign exists 
for the odd-parity case^^^) and a gauge transformation 
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(A(r,a) -^ — A(r,a)). The local recovery of symme- 
try does not occur within the core in case II {Lz — 2). 
Thus, the cancellation of angular momentum due to vor- 
ticity and chirality is crucial for this phenomenon. It 
should be noted that such a restoration of symmetry 
occurs only in the pair-potential; for example, there ex- 
ists non- vanishing circular current within the core. How- 
ever, this pseudo time-reversal-symmetry gives an intu- 
itive interpretation on the impurity effects in vortex core. 
From the symmetry and angular momentum, the chiral 
p-wave vortex core with Lz — can be identified with 
the s-wave superconductors in the absence of the mag- 
netic field. For the latter, non-magnetic impurities are 
known to be harmless (the Anderson theorem'*^' ) . These 
two explain why non-magnetic impurities are harmless in 
the chiral p-wave vortex with L^ = 0. 

§6. Conclusion 

The physics of isolated vortices in pure chiral super- 
conductors is similar to that in pure isotropic s-wave su- 
perconductors. The presence of impurities little affects 
the chiral branch in the chiral superconductors. There- 
fore, it is expected that the KP effect in the chiral su- 
perconductors is much more accessible in experiments 
than that in s-wave superconductors. Further, the local 
recovery of modified time-reversal symmetry around vor- 
tex core is found when the vorticity and chirality have 
the opposite sign. 
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